ELSEVIER 


<D 

CrossMark 


Available  online  at  www.sciencedirect.com 

ScienceDirect 


Proceedings  of  the  Combustion  Institute  35  (2015)  1469-1476 


Proceedings 

of  the 

Combustion 

Institute 


www.elsevier.com/locate/proci 


Turbulent  transport  in  premixed  flames 
approaching  extinction 

K.H.H.  Goh* 1,  P.  Geipel 2 ,  R.P.  Lindstedt* 

Department  of  Mechanical  Engineering,  Imperial  College,  London  SW7  2AZ,  UK 
Available  online  23  June  2014 


Abstract 

The  turbulent  transport  in  premixed  flames  approaching  extinction  has  been  characterised  in  terms  of 
statistical  properties  using  an  opposed  jet  burner  featuring  fractal  grid  generated  turbulence.  The  burner 
was  used  in  a  symmetric  twin  flame  configuration  featuring  pre-vaporised  cyclopentane  and  JP-10  (exo- 
tetrahydrodicyclopentadiene)  at  equivalence  ratios  of  0.75  and  0.85.  The  choice  of  fuels  follows  from  the 
practical  importance  of  JP-10  as  an  aviation  fuel  with  cyclic  C5  compounds  appearing  as  breakdown  prod¬ 
ucts.  The  bulk  velocity  was  set  to  3.0  m/s  resulting  in  a  turbulent  Reynolds  number  of  120.  The  obtained 
data  includes  conditional  velocity  statistics,  flame  curvature  information  and  scalar  fluxes.  The  conversion 
between  conventionally  and  Favre  averaged  statistics  follows  the  Bray-Moss-Libby  theory  and  the  assump¬ 
tion  that  finite  reaction  interface  thickness  effects  can  be  neglected.  The  impact  of  deviations  from  the  latter 
on  statistics  is  explored  and  results  suggest  that  the  effect  is  modest  for  interface  thicknesses  less  than  20%  of 
the  turbulent  flame  brush.  The  experimental  data  obtained  is  sufficient  to  enable  terms  up  to  and  including 
triple  correlations  to  be  evaluated  in  closed  form.  The  work  also  clearly  illustrates  the  rapid  transition  from 
non-gradient  to  gradient  turbulent  transport  as  the  extinction  limit  is  approached. 

©  2014  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Bilger  et  al.  [1]  observe  that  turbulent  premixed 
combustion  is  commonly  agreed  to  be  inherently 
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complex  with  increasingly  wide  application  in 
internal  combustion  engines,  modern  gas  turbines 
for  power  generation  and  jet  engine  afterburners 
[2],  The  increased  use  stems  from  a  desire  to 
reduce  emissions  while  obtaining  higher  efficien¬ 
cies.  Bray  et  al.  [3]  derived  the  arguably  most  suc¬ 
cessful  framework  (BML)  for  analysing  turbulent 
premixed  combustion  on  the  basis  of  a  statistical 
description  and  Bray  [4]  presented  a  recent  review 
of  the  approach.  The  current  contribution  consid¬ 
ers  statistically  stationary  flames  in  the  context  of 
this  framework  through  experimental  measure¬ 
ments  using  a  twin  flame  opposed  jet  geometry 
featuring  fractal  grid  generated  turbulence.  The 
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density  segregation  technique  derived  by  Goh 
et  al.  [5]  is  used  to  detect  flame  iso-contours  and 
to  obtain  conditional  statistics.  Opposed  jet 
flames  have  been  used  extensively  in  the  past 
(e.g.  [6-10])  and  recent  contributions  include 
those  of  Geipel  et  al.  [11],  Goh  et  al.  [12,13]  and 
Coriton  et  al.  [14]. 

The  use  of  fractal  grids  to  generate  multiscale 
turbulence  was  proposed  by  Vassilicos  and  Hunt 
[15]  and  subsequently  investigated  by  Hurst  and 
Vassilicos  [16]  and  Stresing  et  al.  [17],  Fractal  grids 
have  shown  significant  promise  in  generating  ele¬ 
vated  levels  of  turbulence  while  maintaining  flows 
free  of  bulk  instabilities.  The  ideal  grid  for  the 
opposed  jet  configuration  was  identified  by  Geipel 
et  al.  [1 1]  and  found  to  produce  a  100%  increase  in 
Reynolds  stresses  compared  to  traditional  grids 
with  4  mm  holes  as  used  previously  by  Bohm 
et  al.  [18].  Comprehensive  statistical  information 
for  twin  methane  and  propane  flames  has  also 
been  obtained  by  Goh  et  al.  [13]  using  synchron¬ 
ised  velocity  and  scalar  statistics  to  provide  condi¬ 
tional  dissipation  and  flow  field  structure 
information  via  conditional  Proper  Orthogonal 
Decomposition  (CPOD).  A  revised  configuration 
was  used  by  Goh  et  al.  [12]  to  study  the  transition 
from  conventional  premixed  turbulent  JP-10 
flames  to  a  Homogeneous  Charge  Diffusion  Igni¬ 
tion  (HCDI)  [19]  related  flameless  oxidation  mode 
using  single  JP-10  flames.  It  was  observed  that 
conventional  flames  ceased  to  exist  after  transition 
to  flameless  oxidation,  with  reaction  interfaces 
residing  on  the  instantaneous  stagnation  plane, 
similar  to  flames  close  to  extinction  as  observed 
by  Kostiuk  et  al.  [6].  Furthermore,  a  shift  from 
non-gradient  to  gradient  turbulent  transport  was 
observed  for  flames  beyond  the  normal  extinction 
point  [12]  as  the  combustion  mode  shifted  to 
flameless  oxidation.  Under  such  conditions,  reac¬ 
tant  conversion  must  be  supported  by  hot  combus¬ 
tion  products  emerging  from  one  of  the  nozzles. 

By  contrast,  the  current  study  considers 
turbulent  transport  in  flames  approaching  extinc¬ 
tion  in  a  conventional  back-to-back  twin  flame 
configuration  that  can  more  readily  be  related  to 
conventional  premixed  turbulent  flame  theories 
(e.g.  BML).  Density  segregation  methods  e.g. 
[5,12,20],  used  in  the  current  work,  result  in  flame 
isocontours  one  pixel  in  width.  For  flames  with  a 
broader  reaction  interface,  the  probability  of 
detecting  a  burning  mixture  may  become  signifi¬ 
cant.  The  effect  has  been  assessed  experimentally 
through  comparisons  with  OH-PLIF  [12,21]  and 
the  current  work  provides  a  further  theoretical 
analysis  of  the  impact  on  statistics.  The  approach 
is  based  on  a  stochastic  analysis  of  the  impact  of 
the  ratio  of  the  flame  thickness  (<5y)  or,  more  gen¬ 
erally,  the  instantaneous  interface  thickness  (<57) 
to  the  turbulent  reaction  zone  (flame  brush)  thick¬ 
ness  (<5, ).  The  analysis  supports  the  conclusion  [12] 
that  the  impact  on  reaction  progress  variable  (c) 


statistics  is  negligible  under  the  current  experi¬ 
mental  conditions.  The  finding  allows  direct  com¬ 
parisons  of  experimental  results  with  theoretical 
investigations  of  premixed  turbulent  flames  stabi¬ 
lised  in  evolving  multiscale  turbulence.  Finally,  a 
set  of  equations  [22]  that  enable  the  closed  form 
conversion  between  conventionally  and  Favre 
averaged  terms,  up  to  and  including  triple  correla¬ 
tions,  on  the  basis  of  presented  experimental  data 
are  reported. 


2.  Experimental  techniques 

The  turbulent  opposed  jet  configuration  corre¬ 
sponds  to  that  of  Goh  et  al.  [12]  with  the  exception 
that  the  optimal  fractal  grids  identified  by  Geipel 
et  al.  [11]  were  used  in  both  nozzles  in  order  to 
produce  a  twin  flame  configuration.  Impact  plates 
were  placed  upstream  of  the  fractal  grids  in  order 
to  isolate  upstream  conditions  from  grid  generated 
turbulence  [23]  and  hence  to  simplify  the  boundary 
conditions  leading  to  a  reduced  size  of  the  physical 
domain  required  for  computational  studies.  The 
nozzle  exits  were  one  nozzle  diameter  D  (=  30 
mm)  apart.  Twin  premixed  flames  of  cyclopentane 
and  JP-10  were  investigated  using  equivalence 
ratios  (4>)  of  0.75  and  0.85.  Flow  rates  were  such 
that  mixtures  from  each  nozzle  had  a  bulk  velocity 
(Ub)  of  3.0  m/s  (at  298  K).  Cyclopentane  and  JP- 
10  were  heated  to  353  K  and  473  K,  respectively, 
with  nozzle  exit  temperatures  of  333  K  and 
408  K  maintained  in  order  to  prevent  fuel 
recondensation. 

The  flow  control  system  was  the  same  as  used 
by  Goh  et  al.  [12]  with  uncertainties  in  flow  rates 
^0.8%  for  each  fluid.  Flow  rates  of  cyclopentane 
and  JP-10  were  metered  and  vapourised  using 
Bronkhorst  CoriFlow  M53  and  CEM  W-303A 
units  with  the  vapourised  fuel  stream  mixed  with 
air  and  split  equally  using  two  needle  valves.  Mea¬ 
surements  performed  using  GCMS  on  the  reac¬ 
tant  stream  using  an  Agilent  7890A  series  GC 
with  a  5975C  inert  MSD/DS  Turbo  El  Bundle 
equipped  with  60  m  DB-1  column  have  confirmed 
that  no  cracking  of  the  fuel  takes  place.  Coflow 
velocities  were  set  to  0.3  m/s  to  remove  any  large 
scale  bulk  motion.  The  reactant  streams  were 
mixed  further  with  seeded  air  within  heated  hoses, 
temperature  controlled  with  an  accuracy  of  ±1  K, 
for  about  3  m  before  being  introduced  into  the 
nozzles.  A  digital  camera  was  used  to  capture 
the  CH  chemiluminescence  from  both  flames 
and  measurements  performed  only  when  CH 
intensities  were  equal. 

Particle  Image  Velocity  (PIV)  measurements 
were  conducted  in  2D  using  a  120  mJ  Solo-New 
Wave  Nd:YAG  laser.  Both  the  upper  and  lower 
streams  were  seeded  with  aluminium  oxide 
particles  of  about  3  pm  diameter  with  correlation 
between  Mie  scattered  images  calculated  using 
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decreasing  interrogation  window  sizes  from 
128  x  128  via  64  x  64  to  32  x  32  with  a  50%  over¬ 
lap  resulting  in  velocity  vectors  spaced  about 
0.4  mm  apart.  Laser  pulses  were  separated  by 
40  ps  to  minimise  spurious  vectors,  and  1000 
image  pairs  were  obtained  at  around  5  Hz  for 
each  case  (e.g.  [8,11,12]). 


3.  Postprocessing  techniques 

For  every  image  pair  the  density  segregation 
method  of  Goh  et  al.  [5]  was  used  to  detect  flame 
isocontours  in  order  to  synchronise  the  instanta¬ 
neous  reaction  progress  variable  (c)  with  velocity 
vectors.  The  method  has  been  shown  to  be  accu¬ 
rate  to  the  order  of  the  reaction  interface  thickness 
(<5;)  [5]  and  used  by  Goh  et  al.  [12,13,23]  to  obtain 
synchronised  velocity-scalar  statistics.  The  density 
segregation  method  bypasses  the  need  for  an 
additional  laser  for  flame  front  detection  and  over¬ 
comes  problems  associated  with  laser  alignment. 
One  drawback  is  that  detected  isocontours  (reac¬ 
tion  interfaces)  are  one  pixel  wide  (~0.026  mm). 
Reaction  interfaces  are  often  thicker,  as  shown 
by  Goh  et  al.  [12],  and  care  needs  to  be  taken  to 
ensure  applicability  of  the  method.  The  flames 
investigated  here  are  directly  related  to  the  latter 
study  where  the  use  of  OH-PLIF  via  the  method 
of  Kerl  et  al.  [21]  confirmed  the  applicability. 
However,  the  matter  has  also  been  investigated 
theoretically  as  outlined  below. 

Goh  [22]  developed  an  approach  to  ascertain 
the  impact  of  a  finite  mean  reaction  layer  (flame) 
thickness  {3f)  on  measured  flow  field  statistics 
using  a  stochastic  mapping  method.  It  was  shown 
that  the  effects  of  a  finite  interface  thickness  on 
profiles  of  the  mean  progress  variable  were  insig¬ 
nificant  even  when  it  reached  20%  of  the  turbulent 
flame  brush  thickness  (5,).  The  simulated  flame 
brush  was  defined  using  a  Gaussian  pdf  for  flame 
location,  such  that  the  distance  between  the  5th 
and  95th  percentiles  of  its  cumulative  distribution 
function  ( cdf)  was  equal  to  the  measured  value. 
Profiles  of  instantaneous  reaction  sheets  were  also 
assumed  to  be  in  the  shape  of  a  Gaussian  cdf  with 
the  distance  between  the  5th  and  95th  percentiles 
equal  to  the  respective  measured  values  of  Sf  (or 
5/).  A  random  number  generator  [24]  was  used 
to  simulate  the  locations  of  instantaneous  flame 
isocontours  using  the  Gaussian  pdf  of  the  flame 
location  and  instantaneous  profiles  were  mapped 
onto  the  locations  to  obtain  the  first  and  second 
moments  of  the  ‘true’  reaction  progress  variable. 
The  corresponding  profiles  assuming  that 
flame  sheets  were  infinitesimally  thin  (5f  — >  0) 
were  also  obtained.  First  and  second  moments 
of  the  ‘detected’  reaction  progress  variable  were 
obtained  using  a  method  accurate  to  the  order 
of  3f.  The  detected  isocontours  were  assumed  to 
be  uniformly  distributed  between  the  0.5th  and 


99.5th  percentiles  of  the  Gaussian  cdf  for  the  ‘true’ 
instantaneous  flame.  The  infinitesimally  thin  iso¬ 
contours  were  mapped  relative  to  the  true  signal 
in  order  to  obtain  ‘detected’  first  and  second 
moments  of  the  reaction  progress  variable.  The 
same  technique  was  applied  here  to  ascertain  the 
impact  on  the  first  and  second  moments  of  pro¬ 
gress  variable  statistics.  The  former  method  was 
used  to  show  the  relative  effects  on  progress  vari¬ 
able  statistics  if  the  flames  were  assumed  to  be 
infinitesimally  thin,  while  the  latter  method 
included  the  additional  effects  of  uncertainties  in 
the  detection  of  flame  isocontours. 

To  measure  the  instantaneous  flame  interface 
thickness  a  purpose  written  algorithm  was  used 
with  a  Savitzky-Golay  filter  [25]  to  fit  second 
order  polynomials  to  the  detected  flame  isocon¬ 
tours  by  using  10  and  15  neighbouring  points 
(nf„)  in  each  direction,  with  sample  images  shown 
in  Fig.  1 .  The  same  fitting  algorithm  was  also  used 
with  iif„  set  to  5  to  ascertain  the  impact  on  results. 
The  resulting  mean  absolute  deviations  (compar¬ 
ing  nfj,  =  5  and  15)  in  locations  of  fitted  isocon¬ 
tours  (|Arf|  ~  0.003  mm)  and  direction  of  normal 
vectors  to  the  isocontours  (|A0|  ~  2.5°)  showed 
that  the  location  and  orientation  of  the  fitted  iso¬ 
contours  were  relatively  insensitive  to  the  fitting 
parameters.  Subsequently,  bilinear  interpolation 
was  used  to  map  profiles  of  PIV  image  intensities 
normal  to  the  flame  for  all  points  on  the  isocon¬ 
tours  for  each  set  of  1000  image  pairs.  The  profiles 
for  each  image  were  normalised  using  the  mean 
intensity  in  the  reference  windows  used  for  density 
segregation  as  outlined  by  Goh  et  al.  [12],  The 
mean  profiles  of  normalised  PIV  image  intensities 
were  obtained  using  detected  flame  isocontours 
for  each  set  of  measurements  and  the  mean  flame 
interface  thickness  was  obtained  for  each  set  of 
measurements  by  fitting  error  functions  to  the 


Fig.  1.  Sample  PIV  images  for  all  cases,  overlaid  with 
detected  flame  fronts  (lines).  Mapped  local  flame 
normals  (arrows)  in  direction  of  burnt  products  are  also 
shown  for  sample  points  50  pixels  apart  for  clarity.  The 
integral  length  scale  ( Lj )  is  also  shown. 
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profiles.  The  magnitude  of  i5 f  was  defined  as  the 
distance  between  the  5th  and  95th  percentiles  of 
the  error  functions.  The  flame  brush  thickness 
(<5,)  was  also  obtained  in  the  manner  of  Goh 
et  al.  [12],  Values  of  curvature  h  (positive  when 
concave  to  product  stream)  were  obtained  using 
the  above  method. 


4.  Results  and  discussion 

The  obtained  mean  interface  (<5/)  and  flame 
brush  (<5,J  thicknesses  are  presented  in  Table  1. 
The  latter  were  in  the  range  7.5  ^  <5,  (mm)  ^8.5 
and  hot  wire  measurements  in  the  corresponding 
isothermal  flow  produced  integral  length  scales 
(T;)  of  around  4.0  mm  [22],  Accordingly,  the  flame 
brush  is  approximately  twice  the  integral  length 
scale  in  the  reactant  flow.  Values  of  the  flame 
(reaction  layer)  thickness  5f  showed  a  modest 
decrease  as  the  equivalence  ratio  was  increased 
from  0.75  to  0.85,  similar  to  the  measurements 
by  Goh  et  al.  [12],  where  the  thickness  decreased 
from  about  0.78  to  0.58  mm  as  the  equivalence 
ratio  increased  from  0.60  to  0.80.  The  flame  thick¬ 
nesses  were  similar  for  both  JP-10  flames,  proba¬ 
bly  influenced  by  a  slight  asymmetry  at  the 
equivalence  ratio  of  0.75.  The  error  estimation 
technique  outlined  above  was  applied  using  the 
determined  values  of  Sf  and  5,  with  the  detected 
flame  isocontours  randomly  displaced  from  the 
true  isocontours.  The  resulting  normalised  errors 
in  the  turbulent  flame  thickness  (5,),  which  is 
directly  related  to  the  reaction  progress  variable 
(c),  were  found  to  be  around  0.44-0.86% 
(AS,/S,).  The  highest  5f/5t  in  the  current  context 
was  around  12%  for  the  JP-10  flame  with 
<j>  =  0.75.  Even  when  the  detection  of  flame  iso¬ 
contours  had  uncertainties  of  the  order  of  <5/, 
the  expected  deviations  in  c  were  negligible  and 
profiles  of  dd  showed  that  true  peak  values  were 
around  0.22  compared  to  measured  values  of  0.25 
with  minimal  deviations  in  the  profiles  as  shown 
in  Fig.  2.  Hence,  the  impact  of  a  finite  flame  (reac¬ 
tion  layer)  thickness  on  statistics  is  not  significant 
in  the  current  work. 

It  is  not  possible  to  present  all  the  obtained  sta¬ 
tistical  information  and  the  focus  is  placed  on 
data  obtained  along  the  burner  centreline,  where 
V,  Vn  V p,  and  v'c'  are  zero.  Measurements 


Fig.  2.  First  (left)  and  second  (right)  moments  of 
reaction  progress  variable  statistics  obtained  taking  into 
account  finite  thickness  effects,  using  the  method  of  Goh 
[22]  for  measured  values  of  and  <5,  for  JP-10  at 
equivalence  ratio  of  0.75.  Top  row  -  Actual  values 
(solid)  versus  values  obtained  assuming  zero  interface 
(flame)  thickness  — >  0)  (x  with  dashed  lines).  Bottom 
row  -  Actual  values  (solid)  versus  values  obtained  with 
randomly  displaced  infinitesimally  thin  detected  isocon¬ 
tours  (x  with  dashed  lines). 


presented  include  unconditional  and  conditional 
axial  terms  U,  Ur,  Up,  u'u1,  u'rii'r  ,  u'u1  and  u’d , 
the  corresponding  radial  terms  v'v',  v'rv'r  and 
v'pv'p,  as  well  as  the  first  two  moments  of  the  reac¬ 
tion  progress  variable  c  and  dd.  (The  data  enables 
terms  up  to  and  including  triple  correlations  to  be 
evaluated  in  closed  form  [22]).  The  conversion 
between  conventionally  (e.g.  c)  and  Favre  (e.g. 
c)  averaged  statistics,  more  commonly  used  in 
model  formulations,  also  follows.  The  expansion 
ratio  (t  =  —  —  1),  readily  estimated  (e.g.  via  lami- 

Pp 

nar  flame  computation),  is  the  only  unclosed  term. 
Equivalent  equations  have  been  derived  by  Bray 
et  al.  [3]  for  some  of  the  correlations.  The  Favre 
averaged  axial  scalar  flux  can  be  evaluated  by 
re-arranging  Eq.  (2). 


u'd  =  u"d' 


1  +  t 
(1  +  TC)2 


(2) 


Table  1 


Measured  turbulent  flame  brush  (5,)  and  mean  reaction  layer  (flame)  (&f)  thicknesses.  The  uncertainty  in  the  turbulent 
flame  brush  thickness  ( A<5,/ <5 ,)  was  obtained  using  the  stochastic  method  of  Goh  [22],  All  uncertainties  in  presented 
values  illustrate  the  differences  between  upper  and  lower  nozzles. 


Fuel 

<t> 

<5,  (mm) 

<5 f  (mm) 

A<5 ,/d,  (%) 

c5hio 

0.75 

8.13  ±0.47 

0.803  ±0.026 

0.44  ±0.21 

C5H,o 

0.85 

7.65  ±0.19 

0.751  ±0.026 

0.45  ±  0.05 

JP-10 

0.75 

7.39  ±0.06 

0.829  ±0.052 

0.61  ±0.20 

JP-10 

0.85 

7.51  ±0.13 

0.829  ±  0.000 

0.86  ±0.05 
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These  equations  only  hold  for  infinitesimally  thin 
flames.  The  corresponding  triple  correlation  can 
be  obtained  from  Eq.  (4), 

U  =  (\  -  c)Ur  +  cUp  (3) 


u"u"c"  =  c{  1  -  c)  \u2p  -  U;  +  2 U(Ur  -  Up) 

+  (u'pu'p  ~  «)\ 


(4) 


and,  finally,  the  axial  Reynolds  stress  follows  from 
a  re-arrangement  of  Eq.  (5).  The  corresponding 
equations  for  radial  components  follow  by  simi¬ 
larity  and  the  general  expressions  are  listed  in 
Supplemental  Material. 


u'u'  =  u"u"  +  u"u"c" 


(1  +  tc)  l  1  +  ic 


(5) 


Boundary  conditions  in  terms  of  mean  velocity 
components  and  Reynolds  stresses  at  a  distance  of 
2  mm  from  each  nozzle  exit  are  included  in  Sup¬ 
plemental  Material  in  order  to  aid  computational 
studies.  Centreline  profiles  of  mean  axial  velocities 
are  shown  in  Fig.  3,  where  deviations  from 
straight  line  profiles  (isothermal)  increased  with 
equivalence  ratios,  similar  to  previous  findings 
by  Kostiuk  et  al.  [6].  The  corresponding  uncondi¬ 
tional  Reynolds  stresses  are  shown  in  Fig.  4, 
where  peak  radial  components  at  the  nominal 
stagnation  plane  increased  with  equivalence  ratio 
(cf.  Goh  et  al.  [13]).  Reaction  progress  variable 
statistics,  see  Fig.  5,  show  that  flames  move  fur¬ 
ther  from  the  nominal  stagnation  plane  as  the 
equivalence  ratio  is  increased  (cf.  Goh  et  al. 
[12]).  Conditional  velocities  are  shown  in  Fig.  6, 
where  the  reactant  velocities  decreased  slightly 
and  product  velocities  increased  as  the  equiva¬ 
lence  ratio  was  increased.  Conditional  axial  reac¬ 
tant  and  product  Reynolds  stresses  are  shown  in 
Fig.  7.  The  corresponding  profiles  for  the  radial 
components,  see  Fig.  8,  show  that  the  reactant 
Reynolds  stresses  were  almost  constant,  while 
product  stresses  increased  towards  the  nominal 
stagnation  plane,  with  slightly  higher  values  as 
the  equivalence  ratio  increased.  For  JP-10  flames 
at  an  equivalence  ratio  of  0.75,  radial  reactant 


Fig.  3.  Mean  axial  velocity  along  axis  of  burner. 
Cyclopentane  (left)  and  JP-10  (right)  flames  at  equiva¬ 
lence  ratios  of  0.75  (O)  and  0.85  (□). 


Fig.  4.  Mean  axial  (top)  and  radial  (bottom)  Reynolds 
stresses  along  axis  of  burner.  Cyclopentane  (left)  and 
JP-10  (right)  flames  at  equivalence  ratios  of  0.75  (O)  and 
0.85  (□). 


Centreline  [x/D]  Centreline  [x/D] 


Fig.  5.  Mean  first  (top)  and  second  (bottom)  moments 
of  reaction  progress  variable  along  axis  of  burner. 
Cyclopentane  (left)  and  JP-10  (right)  flames  at  equiva¬ 
lence  ratios  of  0.75  (O)  and  0.85  (□). 


Reynolds  stresses  showed  a  peak  before  the  stag¬ 
nation  plane,  possibly  due  to  slightly  different 
equivalence  ratios  between  the  upper  and  lower 
flames.  Results  for  scalar  fluxes,  see  Fig.  9,  show 
large  uncertainties  for  this  case.  The  proximity 
to  global  extinction  (cj>exl  ~  0.7  for  both  fuels)  cou¬ 
pled  with  the  transition  in  flame  propagation 
mode  are  likely  contributors.  Flowever,  it  can  be 
observed  that  as  the  equivalence  ratio  is  increased, 
there  is  a  transition  from  gradient  to  counter  gra¬ 
dient  transport  and  that  the  same  trend  is 
obtained  for  both  fuels.  It  is  imperative  that  com¬ 
putational  models  reflect  such  phenomena  accu¬ 
rately  as  the  propagation  mode  impacts  flame 
dynamics  and  the  stabilisation  point.  Further¬ 
more,  results  show  that  flame  extinction  in  the 
current  twin  flame  configuration  and  the 
transition  to  flameless  oxidation  with  combustion 
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Centreline  [x/D]  Centreline  [x/D] 


Fig.  6.  Mean  axial  reactant  (O  =  Ur)  and  product 
(□  =  Up)  velocities  along  axis  of  burner.  Cyclopentane 
(left)  and  JP-10  (right)  flames  at  equivalence  ratios  of 
0.75  (top)  and  0.85  (bottom). 


Centreline  [x/D]  Centreline  [x/D] 


Fig.  8.  Mean  radial  reactant  (O  =  u'bJ.)  and  product 
(□  =  Vpvp)  Reynolds  stresses  along  axis  of  burner. 
Cyclopentane  (left)  and  JP-10  (right)  flames  at  equiva¬ 
lence  ratios  of  0.75  (top)  and  0.85  (bottom). 


Centreline  [x/D]  Centreline  [x/D] 


Fig.  7.  Mean  axial  reactant  (O  =  u'ru'r)  and  product 
(□  =  u'pu'p)  Reynolds  stresses  along  axis  of  burner. 
Cyclopentane  (left)  and  JP-10  (right)  flames  at  equiva¬ 
lence  ratios  of  0.75  (top)  and  0.85  (bottom). 

stabilised  against  a  hot  product  stream  [12]  are 
both  accompanied  by  a  transition  to  gradient 
transport.  The  mode  of  transport  is  affected  by  a 
number  of  factors,  including  the  expansion  ratio. 
Lindstedt  and  Vaos  [26]  derived  a  relationship 
based  on  the  balance  between  mean  strain  (pro¬ 
duction  term  tensor)  and  mean  pressure  gradient 
effects.  However,  in  the  turbulent  opposed  jet  con¬ 
figuration  the  geometry  imposed  pressure  gradient 
exerts  a  significant  influence  on  the  transition. 
Efforts  have  been  made  to  account  for  such  effects 

[27] ,  though  a  comprehensive  evaluation  of  the 
resulting  model  formulations  present  difficulties 
in  the  current  context,  while  Chen  and  Bilger 

[28]  analysed  the  transition  using  Bunsen  flames. 
Curvature  statistics  for  the  reaction  interfaces 

was  obtained  using  the  Savitzky-Golay  filter  [25] 
based  technique  outlined  above  and  with  three 


Fig.  9.  Mean  axial  scalar  fluxes  (u'c1)  along  axis  of 
burner.  Cyclopentane  (left)  and  JP-10  (right)  flames  at 
equivalence  ratios  of  0.75  (Q)  and  0.85  (□). 


values  of  nfit  (=  5, 10, 15)  to  obtain  fitted  flame 
isocontours.  By  using  a  range  of  values,  it  was 
possible  to  ascertain  the  relative  sensitivity  of 
derived  curvatures  to  the  fitting  parameter  as 
given  in  Table  2.  It  can  be  noted  that  discrepancies 
are  notably  decreased  as  nfil  is  raised  from  5  to  10. 
The  obtained  pdf  of  curvature  is  shown  in  Fig.  10 
and,  again,  it  can  be  observed  that  for  the  two 
higher  values  there  is  good  agreement  apart  from 
at  zero  curvature.  The  latter  discrepancy  is  very 
strong  when  nf„  is  set  to  5  and  a  plausible  explana¬ 
tion  is  that  too  little  local  information  was  avail¬ 
able  to  obtain  accurate  values  for  large  local 
radii  corresponding  to  the  peak  at  zero  curvature. 
Also,  as  the  fundamental  definition  of  curvature 

(h(s)  =  (x'y/f  —  x//y)/(x'2  +/2)3^)  involves  first 
and  second  order  gradients,  the  uncertainties  are 
expected  to  be  large.  Furthermore,  the  ‘optimum’ 
fitting  parameter  is  arguably  dependent  on  the 
local  spectral  content  as  isocontours  are  a  com¬ 
plex  superposition  of  multiple  scales.  Accordingly, 
an  accurate  determination  of  curvature  statistics 
may  require  prior  knowledge  of  the  geometrical 
nature  of  the  isocontours.  Nevertheless,  the  effects 
of  equivalence  ratio  on  curvature  statistics  is 
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Table  2 

Mean  absolute  uncertainties  in  flame  front  location  (|Arf|),  orientation  in  degrees  (|A0|)  and  curvature  (|A/j|). 
Comparison  between  fitted  curves  obtained  using  a  Savitzky-Golay  filter  with  np,  =  5,  10  and  15. 


Fuel 

0 

nflt  Range 

|Arf|  (mm) 

|A0|  (°) 

| AA|  (mm  ') 

C,H10 

0.75 

5/15 

0.00292 

2.42 

0.656 

C,H10 

0.85 

5/15 

0.00305 

2.48 

0.673 

JP-10 

0.75 

5/15 

0.00298 

2.55 

0.658 

JP-10 

0.85 

5/15 

0.00312 

2.56 

0.687 

C,HI0 

0.75 

10/15 

0.00171 

1.18 

0.154 

C5H10 

0.85 

10/15 

0.00179 

1.21 

0.161 

JP-10 

0.75 

10/15 

0.00177 

1.28 

0.157 

JP-10 

0.85 

10/15 

0.00184 

1.25 

0.165 

Fig.  10.  Curvature  pdf  at  resolution  of  0.1  (1/mm)  for 
cyclopentane  (left)  and  JP-10  (right)  flames  at  equiva¬ 
lence  ratios  of  0.75  (top)  and  0.85  (bottom).  Savitzky- 
Golay  filter  with  np,  =  10  (dashed)  and  15  (dotted). 
Empirical  Gaussian  fit  by  Bradley  et  al.  [29]  is  also 
included  for  comparison  (solid). 

modest  even  close  to  the  global  extinction  limit. 
The  latter  is  consistent  with  related  experimental 
studies  for  stable  flames.  Bradley  et  al.  [29]  ana¬ 
lysed  a  wide  range  of  flames  and  showed  the  appli¬ 
cability  of  a  symmetric  Gaussian  form  of  the  pdf 
of  curvature  with  the  peak  value  at  zero  curvature 
linearly  related  to  the  Damkohler  number  as 
Da1/2.  Yuen  and  Giilder  [30]  also  showed  a  sym¬ 
metric  Gaussian  form  of  the  pdf  for  methane 
and  propane  flames,  while  Gashi  et  al.  [31] 
addressed  comparisons  between  2D  slices 
obtained  experimentally  and  DNS  studies.  Sadan- 
andan  et  al.  [32]  showed  that  a  broadly  Gaussian 
pdf  subject  to  local  conditions,  was  obtained  for  a 
range  of  stable  syngas  flames  operating  on  a  swirl 
burner  at  pressures  up  to  20  bar.  Here,  the  sugges¬ 
tion  by  Bradley  et  al.  [29],  see  Supplemental  Mate¬ 
rial,  results  in  a  peak  value  of  0.73  and  a 
somewhat  narrower  distribution  (see  Fig.  10). 
The  current  study  thus  suggests  that  a  (near) 
Gaussian  distribution  extends  to  flames  undergo¬ 
ing  a  flame  propagation  mode  transition  close  to 
global  extinction. 


5.  Conclusion 

The  turbulent  transport  in  premixed  flames  of 
JP-10  and  cyclopentane  approaching  extinction 
has  been  characterised  in  terms  of  statistical 
properties  using  a  twin  flame  opposed  jet  burner 
configuration  featuring  fractal  grid  generated  tur¬ 
bulence.  A  comprehensive  set  of  statistical  infor¬ 
mation  has  been  provided  along  the  axis  of  the 
burner  to  support  comprehensive  comparisons 
with  computational  studies.  The  obtained  data 
includes  conditional  velocity  statistics,  flame  cur¬ 
vature  information  and  scalar  fluxes  and  permits 
the  closed  form  conversion  between  convention¬ 
ally  and  Favre  averaged  statistics  on  the  assump¬ 
tion  that  finite  reaction  interface  thickness  effects 
can  be  neglected.  It  has  been  shown  that  the 
impact  of  the  latter  on  measured  statistics  is  not 
significant  in  the  current  work.  Results  also 
clearly  illustrate  the  rapid  transition  from  non¬ 
gradient  to  gradient  turbulent  transport  as  the 
extinction  limit  is  approached.  The  transition  in 
flame  propagation  mode  is  important  as  it  affects 
flame  dynamics  and  the  flame  stabilisation  point. 
The  results  also  confirm  that  a  (near)  Gaussian 
pdf  of  flame  curvature  prevails  close  to  global 
extinction  and  that  the  change  in  flame  propaga¬ 
tion  mode  is  related  to  the  disappearance  of  the 
conventional  premixed  turbulent  flame  with  a 
transition  to  a  flameless  oxidation  mode  required 
to  support  combustion  in  leaner  mixtures. 
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Appendix  A.  Supplementary  material 

Supplementary  data  associated  with  this  article 
can  be  found,  in  the  online  version,  at  http:// 
dx.doi.org/10.1016/j.proci.2014.06.030. 
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